Scaling and Renormalization Group in Replica Symmetry Breaking space: Evidence 
for a simple analytical solution of the SK model at zero temperature 
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Using numerical self-consistent solutions of a sequence of finite replica symmetry breakings (RSB) 
and Wilson's renormalization group but with the number of RSB-steps playing a role of decimation 
scales, we report evidence for a non-trivial T — > 0-limit of the Parisi order function q{x) for the SK 
spin glass. Supported by scaling in RSB-space, the fixed point order function is conjectured to be 
q*{a) = ^1 er/(|) on < a < oo where f -> a at T = and ^ 1.13 ± 0.01. ^ plays the role of 
a correlation length in a-space. q*{a) may be viewed as the solution of an efi^ective ID field theory. 

PACS numbers: 75.10.Nr,75.30.Cr,71.23.An 
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It is now 30 years since, stimulated by the seminal 
work of Edwards and Anderson Q [23| , one of the authors 
(DS) devised the canonical soluble mean-field spin glass 
model |22| now known as the Sherrington-Kirkpatrick (or 
SK) model. The first published report already exhib- 
ited the 'smoking gun' which a few years later led Parisi 
to his remarkable predictions of a hierarchy of replica 
symmetry breaking (RSB)0| and its associated implica- 
tions for macrostate complexity Q in disordered frus- 
trated systems. These studies further stimulated new 
fields of investigation across physics, computer science, 
biology, econophysics and probability theory . There 
has been great progress since then, in studies of the origi- 
nal model and its many extensions; the concept of RSB is 
now well-accepted for mean field models |2Jl and the be- 
haviour of the Parisi order function is well studied pertur- 
batively near the transition temperature; there exist self- 
consistency equations from which the Parisi function can 
be obtained in principle; and some features of RSB are 
considered known at least semi-quantitatively down to- 
wards zero temperature. However, here we demonstrate 
that there appear to be unappreciated subtleties in the 
limit as zero temperature is approached. In particular, 
we exhibit evidence for a new zero-temperature limit- 
ing Parisi order function with a new 'correlation length' 
in RSB-space and also unconventional behaviour of the 
T = local field distribution in the zero field limit. 

The models we consider are the usual SK one 



H — ^'^Jijcnaj; cr = ±1 
and the Fermionic Ising Spin Glass (FISG) [H [l^ 
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in the half-filling limit, where the spin variables are given 



by 



\a\ and at and a, are fermion creation 



"i"t 



and annihilation operators. |25j . The Jij are quenched 
disordered and frustrated interactions independently dis- 
tributed over all pairs of sites as 



We measure temperature in units of J and for nota- 
tional convenience set the energy scale J — 1. We shall 
consider its properties in replica theory and assume the 
Parisi ansatz as an infinite sequence of replica symmetry 
breaking (RSB) steps for an order function q(x) defined 
on < a; < 1 0, ■ Denoting the steps as Xk it is be- 
lieved [mil that in the hmit T ^ and smaU x < 0{T) 
the steps scale with T and the natural variables are T- 
normalized steps Uk given by 



ak = lirriT^a {J/T)xk{T), 



(4) 



which become selfconsistently distributed on the interval 
< a < oo. We focus on the T = limit and express 
the free energy functional of the Parisi scheme at a-RSB 
in terms of the variational parameter set {qk+i,ak',k — 
1...Q;}; we also have constrained {qi = l,qa+2 = 0} (261] . 
The extremization has been performed for a = 1...5 to 
an accuracy of 4 decimal places. The resulting values of 
qk,ak are shown in Fig^for a = 2- to 5-RSB, together 
with fitting functions of the form 
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P(J) = exp{-NjfJ{2J^))/{^2^^jyN). 
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where fa{a) is chosen as (a^ -f- w{a)) in a 2-parameter 
model and as {a^ + c{a)a + w{a)) in a refined 3-parameter 
model. Both choices derive from the ansatz of a Gaus- 
sian distribution of error- functions which (in contrast to 
a single error function) fit well the numerical data. The 
^(a), w{a) and c(a) are viewed as parameters in an RG 
decimation-type fiow space with the a playing the role 
of the degree of decimation. Fig[21 shows these flows to- 
gether with analytic fit functions which extrapolate as 
a ^ cx) to lima^ooi{oi) ~ 1.14 and lima^oow{a) = 
(and for the 3-parameter case also lim,a^ooc[a) = 0). 
Fig. ^ shows the quality of the fits. Moreover, if the 
indices k are considered as potentially real continuous 
variables, single parameter rescalings Rak and Tak are 
found to map a and q respectively almost onto single 
curves, as Fig|31 shows, with small decreasing deviations 
as a increases. This scaling behaviour in RSB-space of 
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FIG. 1: Order parameter functions of the SK-model at zero 
temperature. Shown as points are the numerical solutions of 
the finite RSB approximations from 2nd (lowest curve) up 
to 5th order, together with fitting functions qmodei{<^) with 
best-fit C(a) and w{a) and also the predicted limiting func- 
tion q*{a) for ^ — — 1.138 is indistinguishable on 
the given scale). Small dots, approximating the smallest 
qk{a),ak{a) by extrapolating the a-dependence beyond the 
calculated 5th order, up to 10th order, confirm the linear rise 
of q*(a). The insert shows the overlap probability distribution 
function 'P*{q) = dC,{q)/dq derived from (?*(a = (,/{! — C))- 



the k-parametrized order parameter function supports a 
single parameter limiting form. Hence we conjecture that 
the infinite RSB limiting form of the T — s- Parisi order 
function is 



*i \ ^ ^ f( ^ 



(6) 



Thus we predict that there is a finite intrinsic correlation 
length ^ in RSB-space. As noted above, our numeri- 
cal results fit £{a = oo) w 1.138. Assuming low T PaT- 
scalingd 0,111 q{x, T) = q{x/T) up to a; = x^axiT) and 
q{T) ~ 1 — 0{T'^) for x > Xmax^ fitting the Gibbs suscep- 
tibility x(0) = limT^o P Jo dx{l - q{x, T)) = da{l - 
q{a)) — ^VtT'? to its expected value x{T = 0) = 1 in 
the a = oo-limit requires ^ = ^/^/tt « 1.128. Apply- 
ing the same scaling assumption to the internal energy 
C7(r) = -i/3 /q dx{l - q^{x,T)) we obtain ^ « 1.1255 
by equating U{Q) and the T — 0-hmit of the free energy 
[27| . The smallness of this length suggests why already 
low orders [a > 2) of RSB produce good results. FigG] 
compares our prediction with the finite- RSB results. 

The overlap probability distribution (OPDF) for fi- 
nite temperatures is given by|0| 'P{q) = dx S{q — 
q{x, T)) = dx/dq, where x{q, T) is the inverse of the orig- 
inal Parisi function q(x,T). The formulation for q*{a) 
which we have determined can be considered valid only 
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FIG. 2: Numerical renormalization group RSB-Row of the 
replica correlation length ^{a) (insert) and of the spread- 
ing parameter w{a) for the 2-parameter model (3-parameter 
model shows an equivalent flow) with analytic fitting curves 
extrapolating to lima^oa^{a) ~ 1.138 and limc,->aaw{a) = 
are included. Note that the evaluation of the points shown 
requires a numerical accuracy far beyond that visible at the 
scale of Fiein 



for low X scaling as T but suffices to give V{q) for q < 
0(1). Rescaling a to a unit interval via a — a/{l+a), 
the contribution of q*{Q to Viq) is T(l - Ci<l)y^'Po{l) 
where V^iq) = < <5(g - 9*(C)) = dC{q)/dq is the 
T — analog OPDF. This function obeys a non-algebraic 
equation yet allows one to extract an analytical result for 
the small-g gap and a g = 1 divergence as 



Po*('? = 0) = 4/^, r*oiq^i) 
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together with a minimum at qmin ~ 0.535. Note, how- 
ever, that the linearity of dependence of q{x) on x only 
as a = x/T is strictly valid only for a < 0(1) and there 
is predicted a plateau of constant q above Xmax of or- 
der 0.5. This yields a delta function in Viq) close to 
q = 1 and approaching q = 1 as T ^ 0. jl^. One 
might note that the original finite-step Parisi equations 
have been considered analytically in the continuum limit 
as the number of steps go to infinity, to yield q{x) as a 
solution of a set of coupled partial differential equations 
with appropriate boundary conditions [MllSj and in prin- 
ciple this could be continued to zero temperature. Un- 
fortunately, however, the practice is more difficult and so 
far not effected analytically. Numerical solution is again 
possible in principle but, since the boundary condition is 
at the transition temperature, it presents accuracy prob- 
lems to reach comparison with the results of finite RSB 
reported here. Hence one is led instead to use the ex- 
trapolation reported above to suggest an underlying 'field 
theory'. To this end we note that q*{a) is a solution of 
the differential equation 
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Qid) = 0, 
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FIG. 5: Rescaled field distributions Pa{h) (top: full scale, 
RSB orders labelled, bottom: gap region zoomed). The Insert 
shows the rescaling parameter Ac as a function of RSB order. 



FIG. 3: Almost perfect mapping of at parameters (top) onto 
a single curve (below) by rescaling with q — + 1 — > Ra{a — k+ 
1) in RSB space; Ra is shown in insert as a function of RSB- 
order. A similar transformation exists for the q parameters. 
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FIG. 4: Quality of fits illustrated by the distance between nu- 
merical data and fit curves for the 2- and 3-parameter models 
(2pm and 3pm, dashed). Labels indicate RSB order. Misfits 
decrease from 0(10"'^) for the 2pm to O(10~*) (3pm) and 
improve with RSB order for both, while the numerical inac- 
curacy of q- and a-parameters is 0(10"''') and O(10~^) resp. 



where a = a/S_. It can also be conveniently expressed as a 
simple generalized equation of motion in a-space through 
a substitution, to yield (dropping the tilde on a) 



iia^dlog{Q{a)) 



da 



with a cubic 'Hamiltonian' 



(9) 



(10) 



Now we turn to the local field distribution 

P{h)=N-^Y.^6{h-Y,J^■J^,)), 



(11) 



again in the limit T ^ 0. The results for P{h) of 0-5 RSB 
show that for any finite RSB there is a gap around h — Q 
which reduces with RSB order but can be made invariant 
by rescaling the field and (inversely) the probability P 
with an a-dependent factor Aq according to 



P„(/l) = A„P„+i(VAa) 



(12) 



with Aq chosen to yield the same gap for each a. This 
results in a set of Pa{h) with good overlap with one an- 
other as h tends towards the gap edge; this is illustrated 
in Fig. |5| Extrapolating to oo-RSB yields the ft, 
fixed point equation 



P*(/l) = AooP*(ftAoo). 



(13) 



The numerical results for the rescaling parameter A show 
an almost linear increase with a, confirming that in the 
oo-RSB limit P*(0) = 0. 

A crude look at the limiting P{h) for h less than of 
order 0.5 shows a linear slope of P{h) of order 0.3, which 
has been suggested as the correct limit for many years 
15, 16, 17, 18, 19] . However a close look at the approach 
to the gap edge suggests that the situation is more sub- 
tle. This is illustrated in Fig. |H| where it is seen that 
the limit of /i ^ has a much smaller slope of approxi- 
mately 0.17 with a slope of order 0.3 taking over only at 
slightly larger h. The figure insert for dP{h) j dh amplifies 
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FIG. 6: Field distribution close to /i = 0; the lines end 

at the gaps in the distributions. The inset shows dP{h)/dh 
with solid lines corresponding to regions of non-zero P{h) and 
dotted lines showing their linear extrapolation to /i = 0. 



this observation as (i) the Hmit of the gap-edge values of 
dP{h)/dh at 0.17 and (ii) a peak in dP(h)/dh around 0.3 
at shghtly higher h 31]. 

In summary, via scahng and a new renormahzation 
group in RSB space, for the first time a surprisingly sim- 
ple analytical model solution q*{a) and a corresponding 
solvable one-dimensional effective field theory are pro- 
posed for spin glasses in the low temperature limit, invit- 
ing the design of approximations which may cope with 
more complicated problems of related but non-solvable 
natures. Similar considerations have raised novel issues 
concerning the pseudogap in the distribution of local 
fields and densities of states. 
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